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Abstract 

Dipole radiation in QCD is derived to the second order in a s . A power-like evolution of 
the spin-singlet heavy quark operators is obtained to the same accuracy. In particular, 
0(a 2 s ) relation between a short- distance low-scale running heavy quark mass and the 
MS mass is given. We discuss the properties of the effective QCD coupling ctf>{E) 
which governs the dipole radiation. This coupling is advantageous for heavy quark 
physics. 
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Theoretical description of heavy flavor decays benefits from a strong hierarchy 
between the mass of the decaying quark and the typical scale of the strong interac- 
tions, m b ( c ) 3> Aq CD . A current level of experimental precision requires an accurate 
treatment of nonperturbative effects even in beauty decays. A consistent genuinely 
QCD-based framework for such a treatment is provided by the heavy quark expansion 
(HQE). The HQE combines the Wilson operator product expansion (OPE) separat- 
ing the physics of low and high momentum scales with the nonrelativistic expansion, 
to treat the nonperturbative effects originating at large distances. It allows a simul- 
taneous precision account for both perturbative and nonperturbative effects. 

An important class of applications of the HQE are semileptonic weak transitions 
between b and c quarks, i.e. the situation when both initial and final quarks are 
heavy. In particular, the most informative predictions can be made for the transition 
amplitudes in the so-called small velocity (SV) limit 0, when the velocities of heavy 
hadrons in initial and final state are small. The heavy quark in this case plays a 
role of a static (or slowly moving) source of the color Coulomb field which affects 
the light degrees of freedom in the hadrons. This physical picture is formalized in 
the Wilsonian field-theory approach by integrating out high-momentum degrees of 
freedom of full QCD and resorting to an effective low-energy theory, where highly 
excited hadronic states are not present. 

A peculiarity of the effective theory for heavy flavor hadrons is its essentially 
Minkowskian nature. In QCD any process in which the velocity of the heavy quark 
Q changes, involves actual gluon radiation with energy and momentum in the whole 
range up to the quark masses. On the other hand, effective theories are called upon to 
eliminate all high-momentum physics. Such a problem does not arise in usual effective 
theories (typically formulated in the Euclidean space) where all short-distance physics 
is virtual. This necessitates a careful control over the radiation effects. 

On the other hand, when the energy loss is small compared to mq, the radiation 
off the heavy colored particles is almost a classical effect, and therefore a universal 
process-independent description is possible. In the limit of small velocity it is a 
familiar dipole radiation. It has, of course, some peculiarities in the non-Abelian 
theories like QCD. 

Even though the dipole radiation has the most obvious manifestation in the SV 
processes, it is relevant in a more general context, for example, for zero- recoil processes 
or inclusive decay widths. The reason is that the OPE ensures that the effects of 
soft physics originating at the momentum scale well below ttiq enter all genuinely 
short-distance observables in a universal way, via its contribution to the local heavy- 
quark operators. The knowledge of the perturbative dipole radiation allows one to 
determine the short-distance evolution of a number of composite operators in the 
effective theory of heavy quarks. 

In the present paper, we derive the non-Abelian dipole radiation to the second 
order in a s and use this result to obtain a power evolution of a number of spin- 
singlet heavy quark operators to order a 2 s . In particular, a gauge-invariant relation 
between a short-distance low-scale running heavy quark mass and the MS mass 
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J71q(j71q) is given, which becomes a practical necessity since the two-loop accuracy in 
the processes with heavy quarks is becoming the state of the art. 

The radiation by the heavy charged particle occurring when its velocity changes 
is a well-known effect from classical electrodynamics. It is most simply obtained by a 
direct computation of the Lienard-Wiechert retarded potentials A^r, t) at r — > oo, 
and has the form 

l_djp_ = a/1 ln i±M_ 2 \l = ?a^ + 

W d(J 7T \ \V\ 1 — \V I y W 3 7T C<J 

where a; is the radiated energy, is intensity and a is the fine structure constant. 
We will be interested only in the dipole term oc v 2 in QCD, and do not consider 
multipole radiation proportional to higher powers of v. 

In quantum electrodynamics the same relation holds, with 1/u dl/du giving the 
probability to radiate soft photon(s) with energy uo. Moreover, there are no higher- 
order corrections in a em , provided u> <C m and the recoil effects of the radiating 
particle m are neglected. Nontrivial corrections are suppressed by powers of u/m. 

A certain type of corrections emerge only due to effects of polarization of the 
quantum vacuum by dipole radiation. It is, therefore, directly related to the running 
of the gauge coupling in the quantum field theory. The radiation of the real (on-shell) 
photons is proportional to the physical value of the fine structure constant a = a em (0). 
However, because of the pair production, the radiated energy is governed by the 
running coupling a(u). For small lu, all vacuum polarization effects are suppressed 
by powers of uj 2 /m 2 where m e is the lightest charged particle. Therefore, the soft 
radiation in QED is defined by a(0) without any correction, in accord with classical 
electrodynamics. This is expected since semiclassical approximation is parametrically 
justified in the limit uo — > 0. 

The dipole radiation in the non-Abelian theory is quite different. Because of 
gluon self-interaction, the effective coupling increases even in the absence of light 
flavors (it also becomes different for different multipoles). Eventually it grows up 
to values of order 1 and the theory enters the strong interaction phase. At such 
energies nonperturbative effects emerge which generate color confinement shaping 
the spectrum of observed hadrons. 

Our main interest lies in the domain where u is large compared to Aqcd (the 
nonperturbative multipole expansion of the color-singlet systems was discussed in 
||). In this regime, perturbative calculations can be performed. Below we introduce 
a few standard notations used for heavy quarks. 

Let us consider a general process of scattering of an external color-singlet weak 
current J with momentum q on a heavy quark Q in the SV kinematics. For simplicity, 
the initial quark is assumed to be at rest. The initial Q and final state Q quarks can 
have arbitrary masses; however, both masses must be large, so that the nonrelativistic 
expansion can be applied. The SV limit v = q/rriQ , \v\ <C 1 is kept by adjusting 
q appropriately. 
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The current J must have a non- vanishing tree level nonrelativistic limit, otherwise 
it can be arbitrary. The most familiar cases are the time-like component of the vector 
current Jq = Q^oQ or the scalar current J5 = QQ. 

The inclusive processes of the scattering on heavy quarks are described by the 
corresponding structure function w(qo, q) which is a sum of all transition probabilities 
induced by J into the final states with momentum q and energy rrtq + q . Using the 
optical theorem, one represents the structure function w as a discontinuity of the 
forward transition amplitude T(q , q) at physical values of q$: 

T(qo, q) = 7^- i d 4 xe-^ x (Q\T J(x) J+(0)|Q> , w(q , q) = 2ImT(g , q) . (2) 

The causality ensures that T(qo, q) is an analytic function of go i n the upper half- 
plane. In the heavy quark limit the spin degrees of freedom become irrelevant, and 
Eqs. (§) assume averaging over spin states (also color, for perturbative calculations, 
etc.). 

The structure functions have a threshold q™ m corresponding to the elastic transi- 
tion: 

I — m r>V 2 

C n = + ^q - rriq - m Q - m Q + — |— . (3) 
The variable u measures the hadron excitation energy in the final state: 

m^v 2 

u = % - C n - q + m Q -mQ- . (4) 



For simplicity, we consider below the case of equal masses, ttIq = tuq, although 
nothing depends on this assumption. 

The nonrelativistic expansion implies that uj rriQ. On the other hand, the 
perturbative treatment is justified at uj ^> Aqcd, and this hierarchy will be assumed 
in what follows. The nonperturbative aspects will be addressed afterwards. Resorting 
to the perturbative calculations, we can use the quark states instead of actual heavy 
hadrons, which was tacitly assumed above. 

The structure function w takes the following form in the heavy quark limit: 



w(u, v) = N5(u) +^^l + (v*) . (5) 



At v = only the elastic peak is present. The excitations, described by the second 
term, appear at the level ~ or v 2 Jl|. The dipole radiation is described by the 

function d(u). 

Motivated by the dipole radiation in QED, we define a coupling a^' by projecting 
Eq. (|5|) on its second term: 

C F ^M = , im l im -^..^M . (6) 
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The denominator in the last ratio is necessary to get rid of the overall normalization of 
the effective nonrelativistic current (renormalized down to the momentum scale near 
uj). In this form the heavy quark limit mq — > oo yields a finite result. It is important 
that the normalization integral includes also the elastic peak, which makes it infrared 
safe at arbitrary v. On the other hand, the exact upper limit does not matter in the 
SV kinematics, since it affects the ratio only by corrections ~ v 2 . 

Therefore, the inelastic structure function in the SV limit determines the effective 
QCD running coupling driving the dipole radiation of gluons. It is a dimensionless 
quantity, and thus is a function of the ratio Aq CD /c<j. The explicit coefficient in Eq. 
is adjusted to make equal to the bare coupling in the Born approximation. 

It is important that the OPE and factorization of the infrared effects^] ensures 
that aW(u) is universal, in the sense that: 

a) it does not depend on the choice of the weak current J probing the heavy quark 
structure, or on the ratio of the quark masses; 

b) as long as the onset of the quark-hadron duality is passed, there is no dependence 
on the particular type of the initial heavy hadron. This property always holds in the 
perturbative analysis. 

The effective coupling obeys the usual renormalization group evolution 

Po = jC A - l n f > ft = 102 - Y nf ■ 

The coefficients fii and higher are scheme-dependent and therefore differ from the 
standard ones obtained for the MS coupling. Note that because of its physical def- 
inition, the evolution of is properly defined when heavy flavor thresholds are 
passed. 

Although the coupling o#") is most natural for heavy quark decays, in practice 
one needs to know its relation to the standard reference coupling a^ s . We calculated 
it to the second order in a 2 using the technique of the Ref. 0: 



7T 7T 



3 J 2 I 6 12. 



a A 2 , /n^3> 



7T 



a MS / e -5/3+lB2^ ^ 2 lqX fn x2 



7T 











12 J 





+ 0{a£) , (8) 



where Ca = N c for the SU(A" C ) gauge group. 

In the latter form we absorbed the effect of running a s into the redefinition of 
the scale at which the MS coupling is evaluated 0. The term ~ Ca represents the 
so-called genuine (non-BLM) second-order effect; it has purely non-Abelian origin. 

1 Physically, it is merely existence of the effective low-energy theory of heavy quarks. 
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We find that its coefficient is not particularly large, —1.685; it is not small either 
as it happens in a number of other observables. Let us note that a similar effective 
coupling a e g governing emission of soft gluons in the ultrarelativistic case v — > 1 has 
been also calculated to two loops ||: 



«eff(/i) . . «f S (/i) 
7T 7T 



6 2 I 12 3, 
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(9) 



its exact definition, however, is less direct and more complicated. It differs from 
in the second-order terms, although the difference is small numerically. 

The BLM-type terms which account for the effects of running of the QCD coupling 
in the leading-order calculations can be easily computed to any order in this case.0 
They contain corrections ~ (j3oa s /n) k where the factor (3q is singled out by the 
dependence on the number of light flavors 0. For example, the BLM estimate of the 
«j term in the last Eq. flD is 

M-"""- (10) 

For book-keeping purposes we also give here the massive quark 0(a s ) 2 contribu- 
tion to a^': with an additional massive quark q the extra term in the last Eq. (J|) is 
given by 



6 V 7T J { 3 v Jx t \ 2/ V 2tJ V t 

where x = 4m^//i 2 and a^ IS (/x) is assumed to be defined for rij + 1 light flavor. 

A construction of the HQE in QCD requires an accurate definition of the basic 
objects of the effective theory, in the first place the heavy quark mass, kinetic oper- 
ator and all other composite operators. Being defined consistently, all they depend 
on normalization point. A normalization-point-independent pole mass of the heavy 
quark, although appearing in the purely perturbative calculations at a given order, 
cannot be completely defined at the level when the nonperturbative effects are ad- 
dressed. By the same token, the parameter A measuring the difference between the 
heavy hadron mass and mq in the heavy quark limit, suffers from a similar uncer- 
tainty. Instead, one must use the short-distance mass mg(/i) and, correspondingly, 
A(/i) = lim^-.oo M Hq - m Q (/i) with Aqcd < i" < rn Q @. The low-scale short 
distance mass can be defined in different ways. We recall, however, that the standard 
MS mass cannot be used for /x rriq [S]. Problems, similar to those in the pole mass, 
also emerge in attempts to define spin-singlet operators without a powerlike mixing 
with lower- dimension operators, including the unit operator. 

2 All necessary expressions are given in |J; c n in Eq. (28) is just C'f • o^t^ip)/^ in the BLM 
approximation. The BLM series has a finite radius of convergence, 4//3o (l + -^i) ^ 2 in terms of 
ap(/i), or 4//?o in terms of af^Oe" 5 / 6 ). 
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The physically appropriate gauge-invariant scheme was suggested in || [H| [| and 
is based on the SV sum rules relating the moments of the SV structure functions to 
the local heavy quark operators. The normalization point is introduced as an upper 
cutoff in the integral over excitation energy. For example, for a heavy hadron Hq 
one defines 



AM 



lim 

mg-»o 



M, 



H, 



m Q ((x) 



lim lim 



Jo 



/40) 



[H Q \Q{iDfQ\H Q ), 
(H Q \QQ\H Q ), 



3 Jq uj 2 w(lu, v) dui 



lim lim 



(12) 



(13) 



etc. The above operator relations for the SV moments of the structure function 
suggest that the perturbative evolution must be driven by one and the same function. 
With the above definitions it is directly a^Qj,). From Eqs. (|1^,|T3|) one finds: 



dA(/i) 
d/i 



4 a {u}) 

f^l s 

o^F 

3 7T 



d/i 2 



a 



(A*) 



7T 



(14) 



where of^(^t) is given in Eq. @. 

We note that the right-hand side of Eq. (|14]) for the kinetic operator was obtained, 
with the a 2 accuracy, in JIT] considering the zero-recoil processes, where the dipole 
radiation is absent and the inelastic structure functions appear only in the order 
l/m/Q. Once again, the OPE ensures that the result is process-independent, which is 
demonstrated here at the quite nontrivial level of genuine two-loop corrections.^ The 
SV limit, however, is technically important since the separation between high- and 
low-energy hadronic states by means of the cutoff in u is velocity-independent only 
through order v 2 . 

The similar /z-dependence holds for the renormalized slope of the Isgur-Wise func- 
tion when it is defined in the same physical way ||, based on the Bjorken sum rule: 



d/i 



4 2a^(fi) 
3 3tt 



(15) 



Due to the short- distance effects, the observable transition form factors do not have 
a literal heavy quark limit at v ^ but vanish, and require factoring out the per- 
turbative suppression due to gluon radiation, to yield the effective /i-dependent IW 
function £(t> 2 ; /i). The renormalization factors generally differ depending on the way 
the /x-cutoff is introduced in the perturbative factors. In analogy with Eq. ([I]), it is 
convenient to define £(t> 2 ; jj) for small v as 



Z(v 2 - y) 



lim 

ruQ — >oo 



F(v) 



2^ Jo v) du 



1/2 



f? 2 M = "2 



dv 2 



\v=Q 



(16) 



5 To first order in a, it was checked in 



and, for the BLM terms, in fl6[ to all orders. 
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where F is the form factor, for example, of the vector current Q^qQ and w is the cor- 
responding structure function. This is the scheme suggested in || for which Eq. ( |l5f) 
holds. Various definitions using the dimensional subtraction schemes differ from this 
one; they do not have such a direct physical meaning, however. It is interesting to 
note that the perturbative effects in g 2 (fi) can be viewed as a universal /^-dependent 
renormalization of the tree-level constant accompanying it (—1/4 for scalar mesons 
or for baryons). 

The RG equations ([14]) are directly proportional to a.f\ Hence, the application of 
the OPE shows that the first-order running is not renormalized in the Abelian theory 
(without light fermions) to all orders in the coupling ||, since in such a theory a^' 
identically coincides with a(0), as discussed above. 



According to fi~0|l , one can calculate perturbatively a physical (gauge-invariant) 
low-scale running mass mg(/i) with fi <C niq by subtracting, order by order, the 
infrared part given by the SV sum rules, from the pole mass mg le : 



m 



polo 



pert 



[Aqcd(^)] 



pert 



2™q(^) 



v 2 M 



pert 



This mass determines the kinetic energy term in the renormalized heavy quark Hamil- 
tonian when it is expanded in the velocity of the heavy quark: 



H Q = QA Q + 




;i7) 



The two-loop relation between this mq(n) and the MS mass is 



+ 



a s [m) 



7T 



_ J 4 a s (m) 
mil + 

I 3 7T 



2 r 



K 



3 m 9 




zz ] + i£ [ hx 2 + — 

m 2 \ 3 3 



where 



K 




X 1 



665 
144 



m 3 



71 

" 18 

Xo 



21n2 



19 

Y 



m 



, 2a 13 
In 
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(19) 



(20) 



and m = m{m) is the MS mass normalized at the scale m. We neglected small terms 
~ yU 3 /mg which can be incorporated in the same way, if necessary. We used here the 
0{a 2 s ) relation between the pole mass and the MS mass ||12|| . 

The BLM part of the a z s terms in mg(/i)/m in Eq. (|18D is given by 



m 




{ 2353 13 2 7 >/ N 16 a 

i + — 7T 2 + -C 3 — 

[2592 36 6 s v ' 9 m 



67 



7T~ 



xt + 

1 36 6 
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2^_ 

3 m 2 



xi + 



10 _ TT^ 

T ~ If 




(21) 



where the expressions of [[/J are used. 

A word of clarification is in order regarding the definition of a^K We did not 
specify in Eq. (|J) which heavy quark masses should be used to define the reference 
point for uo. Clearly, they must be low-scale masses, determined in one and the 
same way for both initial and final state quarks. The dependence on the choice then 
disappears in the difference (nig — tuq) in the heavy quark limit, but still persists in 
the kinetic energy of the recoiling quark given by the last term in (^). 

Since a possible difference in the definition of uo by itself is quadratic in velocity, 
the choice of the normalization point for the quark masses does not matter herefj This 
is, again, the reason why the analysis becomes so much simpler in the SV kinematics. 

We conclude with addressing the question of nonperturbative effects in the dipole 
radiation. When the effects suppressed by powers of Aq C d/^ are not neglected, the 
nonperturbative physics starts to play a role. In particular, the actual structure 
functions in Eqs. (Q) refer to a particular heavy hadron (B meson, A& etc.) but 
not an isolated heavy quark. Likewise the complete effective coupling ac£ > ' defined 
in this way, as any other effective charge [f[4|j , acquires certain power-suppressed 



nonperturbative contribution. For example, this coupling would differ in B and A& 
even in the heavy quark limit. 

One can quantify the nonperturbative effects in a^'(u) by constructing its 1/uo 
expansion. f] We will show now that the methods used in the OPE are applicable here. 

In general, the transition amplitude T(oo, v) has a number of cuts. The physical 
'decay' cut, which generates w(u,v), starts at uo ~ 0. The other cuts go away to 
infinity when ttIqq — > oo ||. Therefore, a^\w) is completely determined by the 
analytic continuation of the nonrelativistic part of T(uo,v) to positive uo. In the 
dispersion relation 

m/ 1 f w(uj'v) , , If discT(uo'v) , . 
T(u, v) = — / - V . duo' + — / — K —^du' (22) 

Ju'>Q UO — UO — 16 2m Jother cuts UO 1 — UO — IE 



2tt 



the contribution of the additional cuts vanishes in the heavy quark limit. At the same 
time, the large-cu behavior of T(uo) in the complex plane and, therefore, the asymp- 
totics of a^\uo)/uo is determined by the standard OPE expansion which expresses 
the contribution of the soft physics in terms of the expectation values of local heavy 
quark operators QOkQ of dimension k + 3: 

So^H _ ^ 1 (H Q \QO k Q\H Q ) 

uo ~\ Ck 2M HQ uo^ ■ {26) 



4 One should, actually, take the hadron masses in the kinematics. 

5 We address here only the effects of the infrared origin and discard a more theoretical problem 
of ultraviolet renormalons ||l5| or effects of small-size instantons. Thus we can merely consider, for 
example, the difference between the coupling from heavy mesons and heavy baryons. 
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In this respect, the case of is similar to the effective couplings defined through 
the Adler D-function a^^Q 2 ) (with D(Q 2 ) = 1 + a[ D \Q 2 ) /n) or a similar ai R \s) 
effective coupling which measures the strong-interaction corrections in the absorptive 
part of the correlator of the vector currents of massless quarks in the Minkowski 
domain. In that cases the OPE predicts that the nonperturbative effects associated 
with the vacuum expectation value of the operator G 2 af5) are suppressed by a fourth 
power of energy. 

The first nontrivial heavy quark operator is the kinetic operator Q(iD ) 2 Q, which, 
by dimensional arguments, would lead to 1/u 2 effects. However, it does not contribute 
to ay' . Indeed, we can consider the initial heavy quark moving with the spacelike 
momentum \p\ <C uj. The structure function Wp(uj, v) for such a state \Qp) is simply 
related to the one at rest, wq(lj, v): 

Wp(u,v) = w (uj-vp, v) +0(l/m Q , w 5 ) . 

If we consider a superposition of such states / d 3 pQ(p 2 ) \Qp) which is invariant under 
rotations, its structure function differs from wq only in terms ~ p 2 v A which are 
negligible in the SV limit. On the contrary, the kinetic expectation value for such a 
state is non-zero and equals to (p 2 )- This proves that the corresponding coefficient 
function C2 must vanish. 

With the heavy quark spin decoupled, only spin-singlet operators can appear in 
the expansion Eq. fl23|). There is one spin-singlet, so-called Darwin operator Ob = 
—iQDkD Q DkQ which yields 1 /^-suppressed terms.Q We do not see any reason why 
such terms cannot appear. Moreover, using the sum rule for the Darwin term || 
TBI, O, one finds a unique contribution associated with its logarithmic anomalous 



dimension: 

5oa^{u) 3n%(u)pl(u) ( a s {p) \ WA 3 lD p D {p) /0.55GeV\ 
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) , (24) 



a^\uj) 8ai w) (w) ^ 3 " {"si")] ^ 3 
where 

- 1 \ s 1 \ do.. A , ocf^ in) 13 

IdWPUN = PfyPDW > lD\p) ^ Id 2?r , Id = -y , 

PM = (Bl °°l B) » - 2 ^flM B - 0.1 GeV 3 at a s (p) K 1 . 

We used here the factorization estimate of p D for pseudoscalar (vector) mesons [|16|1 ; 
the anomalous dimensions of the D = 6 operators were calculated in ||18|| . In any 
case, the nonperturbative effects are expected to die out quickly with energy: 



^af)( W ) ~ f^)*, k>3. 



6 At the level of nonperturbative effects, for heavy quark states with spin of light degrees of 
freedom j > 1 there can be the second Lorentz structure in d(u>) ~ (qj) if averaging over these 
spin states is not performed. Then an independent D = 6 operator with two spacelike derivatives 
can appear. 
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Even though is one of many possible effective couplings in QCD, it is most 
appropriate for heavy flavor electroweak transitions. Moreover, it has an additional 
advantage compared to a few alternative ones considered in the literature, like 
or ay 1 '. The latter determine the deviation of the corresponding observables, D(Q 2 ) 
or R(s), from their tree-level value 1 - whereas a^' directly measures the strength of 
the interaction. This is an obvious advantage for its accurate determination. In this 
respect, the two-loop relation for derived in the present letter, has its counterpart 
in three-loop calculations of D-function or R(e + e~ — > hadr). 

From this perspective, is closer to the so-called V^-scheme coupling re- 
lated to the heavy quark potential V(R) |T^, f|, |(J. The effective coupling 



is physically appropriate in quarkonia and near-threshold perturbative calculations. 
However, its exact definition in higher orders is not evident at the moment. As a 
matter of fact, it requires knowledge of the heavy quark potential at all distances R 
including R — > oo, even when ct^\q 2 ) is evaluated in the perturbative domain. In 
particular, appearance of terms a k s In a s with k > 4 in its relation to the [ITS 



may indicate that this is not a genuinely short- distance quantity. In contrast, crp> 
is completely defined, not only in the perturbation theory but even beyond it. This 
clearly represents a certain theoretical advantage. 

One aspect of ay' must be noted: it is defined in Minkowski space and in this 
respect is similar to ay® but not to a^ s D \ Since, as was mentioned, in this problem the 
Euclidean OPE is still applicable, it is not too important. Nevertheless, the so-called 
"exponential" terms in related to the asymptotic nature of the power expansion 
of the OPE, which are truly exponentially suppressed in the Euclidean domain, may 
yield only oscillating power-suppressed component in a[ R \ and, similarly in ay\ The 
recent theoretical discussion can be found in [[H]. The coupling ay' may also show 
some (smeared) resonance structure at the heavy flavor thresholds. 

The Euclidean counterpart of ay' can also be defined in a straightforward way. 
One can either use a Borel transform image of the structure function a^(u)/uj, or 
usual dispersion integral. These technical details will be given elsewhere. 

The MS coupling a^ s proved to be indispensable for complicated multi-loop cal- 
culations. On the other hand, it is rather unphysical in some respects and its use as 
an expansion parameter sometimes obscures the perturbative expansion | 22|j . At the 
moment the complete 0{a 2 s ) relation between a^ s and affl is enough. The order- 
ed calculation seems feasible as well.Q It is natural to think that using ay 3 ' or its 
Euclidean counterpart is advantageous in heavy quark decays, and may improve the 
accuracy of the perturbative estimates. 

To summarize, we derived the non-Abelian dipole radiation by nonrelativistic 
color particle to second order in the strong coupling. By virtue of the SV sum rules 
this determines the two-loop power mixing of the number of spin-singlet effective local 
heavy quark operators. We give the two-loop relation between the gauge-invariant 



7 It is interesting to attempt calculating the perturbative corrections to ai by considering dipole 
radiation in the nonlinear non-Abelian classical colordynamics. 
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short- distance low-scale heavy quark mass mg(/i) suitable for the Wilson OPE and 
the value of m^mg) in the MS scheme. We discuss the properties of the effective 
dipole radiation coupling a^' which is useful for heavy quark physics. 
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